Anomalous Polymer Dynamics Is Non-Markovian: Memory Effects and The 
Generalized Langevin Equation Formulation 
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Any first course on polymer physics teaciies tiiat tlie dynamics of a tagged monomer of a polymer 
is anomalously subdiffusive, i.e., the mean-square displacement of a tagged monomer increases as 
t" for some q < 1 until the terminal relaxation time r of the polymer. Beyond time r the motion 
of the tagged monomer becomes diffusive. Classical examples of anomalous dynamics in polymer 
physics are single polymeric systems, such as phantom Rouse, self-avoiding Rouse, self-avoiding 
Zimm, reptation, translocation through a narrow pore in a membrane, and many-polymeric systems 
such as polymer melts. In this pedagogical paper I report that all these instances of anomalous 
dynamics in polymeric systems are robustly characterized by power-law memory kernels within 
a unified Generalized Langevin Equation (GLE) scheme, and therefore, are non-Markovian. The 
exponents of the power-law memory kernels are related to the relaxation response of the polymers to 
local strains, and are derived from the equilibrium statistical physics of polymers. The anomalous 
dynamics of a tagged monomer of a polymer in these systems is then reproduced from the power- 
law memory kernels of the GLE via the fluctuation-dissipation theorem (FDT). Using this GLE 
formulation I further show that the characteristics of the drifts caused by a (weak) applied field on 
these polymeric systems are also obtained from the corresponding memory kernels. 

PACS numbers: 05.40.-a, 02.50.Ey, 36.20.-r, 82.35.Lr 



I. INTRODUCTION 

One of the most common phrases of wisdom in polymer 
physics is that a polymer, in its terminal relaxation time 
T, displaces itself by its own size in physical space [J 
Being constructed from monomers connected in series, a 
polymer has a wide range of length-scale dependent re- 
laxation times. The longest one of them is the terminal 
relaxation time r, manifested in the decay of the poly- 
mer's end-to-end vector correlation function. [E.g., for a 
phantom Rouse polymer of length N , the relaxation time 
corresponding to length scale N/p scales ^ {N /pf' with 
p a positive integer, implying that r ~ N"^ for a phantom 
Rouse polymer] [2]. The terminal relaxation time r of a 
polymer scales with its length as a power-law t ^ for 
some exponent k, while its own size in space scales as iV" 
for some exponent v. If the polymer is to displace itself 
by its own size in time t, in simplest of cases the mean- 
square displacement of a tagged monomer of a polymer 
must behave ~ t^'^/'^ , and since the quantity 2i//k is not 
necessarily unity, the dynamics of a tagged monomer in a 
polymer must be anomalous till the terminal relaxation 
time. If asked to provide canonical examples of such 
anomalous behavior, a polymer physicist would almost 
certainly cite single polymer dynamics, such as phantom 
Rouse (y = 1/2, k = 2), self-avoiding Rouse \v sa 0.588 in 
three dimensions (3D), v — i/A in two dimensions (2D) 
and K = 1 -|- 2u\ and self-avoiding Zimm [v w 0.588 in 
3D, = 3/4 in 2D and k — iv) polymers. (The correct 
single polymer dynamics in a fluid was first presented 
by Zimm In this paper I refer to such polymers, 

for which the monomers interact with each other via hy- 
drodynamic interactions, as Zimm polymers. Few years 
earlier than Zimm, Rouse [3| put forward a model for 



single polymer dynamics that neglect the hydrodynamic 
interactions between the monomers. Although incorrect, 
the corresponding (Rouse) polymer dynamics, for poly- 
mers that are self-avoiding, as well as the phantom ones, 
i.e., polymers that can intersect themselves, resides at 
the heart of polymer physics Q — and widely used till 
today — largely due to its simplicity.) Given such an 
abundance of anomalous dynamics in polymeric systems, 
a natural question is whether it is possible to provide a 
generic stochastic foundation for it. 

In a recent Letter (5(, hereafter referred to as Ref. I, 
I reported that the anomalous dynamical behavior for 
phantom Rouse, self-avoiding Rouse and Zimm poly- 
mers, and that of polymer translocation through a nar- 
row pore in a membrane can be theoretically formulated 
via the following Generalized Langevin Equation (GLE), 
wherein the velocity v{t) of a tagged monomer and the 
force (pit) it experiences, are related to each other via 



(t>{t} = I dt'fi{t - t')v{t') + g{t). 
Jo 



(1) 



In Eq. ID) fi{t) is the memory kernel, and the stochastic 
noise term g{t) satisfies the condition that {g{t))o — 0, 
wliile the fluctuation-dissipation theorem (FDT) {g{t) ■ 
g{t'))o = {${t) ■ (f{t'))^=o = 3kBTfi{t - t') in 3D. Here 
ks is the Boltzmann constant, T is the temperature, and 
(. . .)o denotes an average over the stochastic noise real- 
izations, including an average over equilibrium configu- 
rations of the polymers at < = 0. Equation ([1]) can be 
inverted to write 



v{t) = - / dt'a(t 
Jo 



t')${t') + hit), 



(2) 
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polymeric system 


scaling of r 


mean square displacement 


flit) 


phantom Rouse 


~ N'^ 


~ t^''^ till T and ~ t thereafter 


~t-^/2exp(-t/r) 


self-avoiding Rouse 




~ till r and ~ f thereafter 


~ t-2-/(i+2-) exp(-t/r) 


self-avoiding Zimm 




~ t^''^ till r and ~ t thereafter 


~t-2/^exp(-t/r) 


phantom Rouse 
polymer translocation 


~ 

(r not translocation time) 


~ t 


'--^ t"""" exp(— t/r) 
translocation time ~ A*'^ 


self-avoiding Rouse 
polymer translocation 


(r not translocation time) 


~ t(n-)/{i+2-) till r 
and ~ t thereafter 


~ t-d+^j/d+a-) exp(-t/r) 
translocation time ~ N'^^" 


self-avoiding Zimm 
polymer translocation 


(r not translocation time) 


~ t(i+-)/{3-) till r 
and ~ t thereafter 


translocation time ~ A*'^^^" 


reptation; repton model 
(curvilinear co-ordinate) 




~ t^''^ till r and ~ t thereafter 




melt (reptation theory) 


~ Ar3 


~ t^/'' between and ~ iV^ (*) 
~ t^/^ between ~ A^^ and ~ r 
and ~ t thereafter 


~ t^^/* between Te and ~ A'"^ 
(as much as data can resolve) 



TABLE I; Summary list of the polymeric systems (polymers of length A') with anomalous dynamics, and the corresponding 
memory kernel fi{t), for which the GLE formulation is reported in this paper. (*) Te is the so-called entanglement time for 
polymer melts in the reptation theory Q. Note that if fi{t) ~ then the anomalous dynamics exponent is also a. 



where jl[s)a{s) = 1 in the Laplace space, {h{t))o ~ 0, 
with the corresponding FDT l$,{t) ■ h{t'))o = {v{t) ■ 
1^=0 = iksTait - t'). In Rcf. I, I argued that 
on the one hand /i(t) is the mean relaxation response of 
the polymers to local strains, and can be derived from 
the equilibrium statistical physics of polymers; and on 
the other, a{t) characterizes the anomalous dynamics via 
the FDT: as the mean-square displacement of a tagged 
monomer is obtained by integrating {v{t) ■v{t')) twice 
in time. An important property of the anomalous dy- 
namics that transpires through this exercise is that if 
fi{t) ^ t^" for some a, then the anomalous dynamics is 
also a. In other words, according to Ref. I, the anoma- 
lous dynamics for phantom Rouse, self-avoiding Rouse 
and Zimm polymers, and that of polymer translocation 
through a narrow pore in a membrane are connected to 
the mean relaxation response of the polymers to local 
strains. 

The purpose of this paper is to provide a pedagogi- 
cal account of the unified GLE foundation (Il][2]) for the 
classical examples of anomalous polymer dynamics, in- 
cluding the ones considered in Ref. I. The precise issues 
covered in this paper are the following. 

(i) I provide an elaborate derivation of Eq. (IT][2]) for 
phantom Rouse polymers, which is the only fully 
analytically tractable polymer dynamical system. 

(ii) I substantially supplement the scaling results of 
Ref. I with simulation data, and extend the GLE 
formulation to the cases of a single reptating poly- 
mer, and polymer melts. 

(iii) Given the GLE formulation I show that the char- 
acteristics of the drifts caused by a (weak) ap- 
plied field on these polymeric systems, too, are ob- 



tained from the corresponding memory kernels: if 
fi{t) ^ t~°' for some a, then the anomalous dynam- 
ics is also a. This could be thought of as the analog 
of the Nernst- Einstein relation [6]. 

(iv) Since the GLE is a description of the trajecto- 
ries in phase space, this paper explicitly brings to 
light the non-Markovian character of the anoma- 
lous dynamics for polymeric systems. A natural 
question that then arises is "Is it possible to pro- 
vide a probabilistic description of these trajecto- 
ries in phase space?" While in this paper I leave 
this question for further research, I note that in 
limited context — for the anomalous dynamics of 
polymer translocation through a narrow pore in 
a membrane — fractional Fokker-Planck equation 
(fFPE) has been put forward in recent times f^. In 
this equation, an extension of the standard Fokker- 
Planck equation, anomalous dynamics is a conse- 
quence of introducing power-law waiting times be- 
fore each jump of the concerned particle, as the 
jump length and the waiting times for any jump 
is obtained from fixed probability distributions, in- 
dependently of their values at previous jumps. It 
is worth emphasizing (as it emerges from the GLE 
formulation elaborated in this paper) that there is 
no power-law waiting time for the movements of the 
concerned monomer; instead if it makes a move at 
any time, there is an enhanced chance to undo this 
move in subsequent times; this is where anomalous 
dynamics in polymeric systems stem from. This 
establishes that the fFPE is unsuitable to describe 
anomalous polymer dynamics (polymer transloca- 
tion included) that are considered in this paper. 

The paper is organized in five sections. Section |TT] is 
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devoted to single polymeric systems. Therein I provide 
(a) an elaborate derivation of Eq. for phantom 

Rouse polymers, (b) the characterization of iJ,{t) as the 
mean relaxation response of the polymers to local strains, 
and (c) the GLE formulation for self-avoiding Rouse, self- 
avoiding Zimm, polymer reptation, and polymer translo- 
cation through a narrow pore in a membrane. In Sec. IIIII 
I provide the GLE formulation for the anomalous dynam- 
ics of a polymer melt, a many-polymeric system. In Sec. 
IIVI I show that the characteristics of the drifts caused 
by a (weak) applied field on these polymeric systems are 
obtained from the corresponding memory kernels. This 
could be thought of as an analogue of the Nernst- Einstein 
relation. Section |V] is devoted to a critique of recent at- 
tempts to model by the fFPE the anomalous dynamics of 
polymer translocation through a narrow pore in a mem- 
brane. The paper is finally concluded in Sec. IVII 

But before I begin, in Table |l] I provide a summary list 
of the polymeric systems with anomalous dynamics and 
the corresponding memory kernel ii(t) for which the GLE 
formulation is argued in this paper. 



II. THE GLE FORMULATION FOR 
ANOMALOUS DYNAMICS IN SINGLE 
POLYMERIC SYSTEMS 

The GLE formulation for phantom Rouse 
polymers 



1. 



Derivation of Eq. {IJ) from the Rouse equation 



Consider a phantom Rouse polymer of length N. It is 
described by the Rouse equation; in continuum represen- 
tation it reads [3| 



7- 



drn{t) , d^rn{t) 



dt 



fn{t), 



(3) 



where r„(t) is the location of the n-th monomer at time 
t, 7 is the damping coefficient of the surrounding fluid, 
and k is the spring constant for the springs connect- 
ing the consecutive monomers. The stochastic force 
fm.{t) satisfies the conditions {fn{t)) = and the FDT 
{fn,a{t)fnx{t')) = 2jkBTSim - Ji)5ax5{t - t'), for a, A = 
{x,y,z). Equation ^ is supplemented by the "open" 
boundary conditions that the chain tension of the poly- 
mer at the free ends must vanish; i.e., {dfn/dn)\n=o — 
{drn/dn)\n=N = 0. 

Since the Rouse equation is linear in r„(t), it can be 
solved to obtain all correlation functions using the mode 
expansion technique Two noteworthy results borne 
out of this exercise are: (a) the terminal relaxation time 
r = ^N^ /(n'^k), and (b) the MSD of the middle monomer 
increases as t^^^ until time r, and only after that time 
the motion of the middle monomer becomes diffusive, 
with the diffusion coefficient scaling ^ It is this 

subdiffusive motion of the middle monomer that I obtain 
from the GLE (IT][2]), which in turn I derive exactly from 



Eq. ([2]). For this problem 4>{t) given by the net force 
it experiences (forces exerted on the middle monomer by 
the ones it is connected to, or equivalently, the chain 
tensions at the middle monomer), i.e., 



(pit) = k 



drnit) 



dn 



drn{t) 



n={N/2) + 



On 



n={N/2)- 



(4) 



Note that using Eq. ^ the Rouse equation for the mid- 
dle monomer can be written as 



(5) 



Given Eqs. ([T][S]); an issue that may come to mind is 
the following. While the GLE formulation ([T][2) has a dis- 
tinct memory kernel — I have shown in Ref. I that ^,{t) ^ 
t^^/^ exp(— t/r), and similarly a(t) ^ i"'^/^ exp(— i/r) for 
phantom Rouse polymers — Eqs. ^SS^ do not: e.g., 
in Eq. ^ the velocity of the middle monomer at any 
time is simply proportional to the forces at the same 
time exerted by the monomers it is connected to (plus 
a (5-correlated random noise). From this it may appear 
that the GLE formulation (IT][2]) is in contradiction with 
Eqs. (I3][S]). In truth however, there is no contradiction 
among Eqs. ([T][5]). First of all, the fact that there is no 
explicit memory kernel in the Rouse equation does not 
mean that a phantom polymer does not have memory. 
As noted in the very first paragraph of the introduction, 
the relaxation time corresponding to a length scale N/p 
scales ~ {N/pY with p an integer (2l| shows that a phan- 
tom polymer does have memory; in fact, below I will 
show that the memory kernel /i(t) is indeed built from 
the long relaxation times for fiuctuations at long length 
scales of the polymer. [In this context, it is worth recall- 
ing the work by Zwanzig [s] , wherein the memory of the 
concerned particle arises due to its coupling to a thermal 
bath of harmonic oscillators, effected via a Hamiltonian. 
For the present case, the fluctuation modes effectively 
play the role of a thermal bath that is coupled to the 
motion of a monomer due to chain connectivity, giving 
rise to the memory kernel /i(i).] Secondly, while Eqs. ([T]) 
and being inverse Laplace transform of each other, 
describe exactly one relation between v{t) and 4>{t), Eq. 

provides the second relation between the two quanti- 
ties, so that out of Eqs. (IT][5]) and (O one can form closed 
equations for both v{t) and (t){t). Indeed, in Sec. Ill A 21 
I use the closed form equation for (j){t) to obtain the ve- 
locity autocorrelation function (of course, it is the same 
velocity autocorrelation function as obtained in Ref. I, 
and as summarized in the introduction of this paper). 

Without further ado, below I derive provide the GLE 
([T|) for phantom polymers from the Rouse equation ([3]). 

Three ingredients are necessary to derive Eq. ([T]). The 
first one of them is the dynamics for both halves of the 
polymer when the middle monomer held fixed at, say, TZ. 
Being a phantom polymer, when the middle monomer 
held fixed, the two halves of the polymer evolve indepen- 
dently of each other. With r^(i) = fvi(i) — 7?, the corre- 
sponding equations for each half of the phantom polymer 
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satisfy the (Rouse) equation 



df 



1- 



dt 



= k 



fn{t), 



(6) 



for n 7^ 0. Equation ([6|) is supplemented by tlie 
"open" boundary conditions that the chain tension of 
the phantom polymer at the free ends must vanish; i.e., 
{dfn/dn)\n=o = {dfn/dn)\n=N = (additionally r q = 
by definition) . In order to analyze the motion of the two 
halves of the phantom polymer I define 



Y^'\t) = ^l! dn sin "('^+^)% -;^^^,(t), (7) 



N 



1 r" 

and Y^''\t) = -— / dn su 



7r(2p+ l)n 
N 



i+Ar/2 



it) (8) 



for p = 0,1, . . ., for the right and the left half, such that 

l)n 



p 



N 



■Sin) 



y(')(i)su 



7r(2p+ l)n 
N 



e(-n). (9) 



Note that Eq. ([9]) is consistent with the open boundary 
conditions, as well as the fact that the middle monomer 
remains fixed. The independent evolution of each half 
satisfies the Langevin equation (LE) 



dY 

lp-gf^~k,Ypit)+fp{t), 



(10) 



for y„ 



]Y, 



{1} vir)- 



, 7p = 2N-f, kp = 27r2fc(2p+ l)ViV, 

(/;) = and {fpa{t)fgx{t')) = lpkBT5{t - t')dp,d^x. 
The second ingredient is the equilibrium averages of Yp 

andY^. First, {Y^'^\t)) = {Y}^\t)) = by isotropy. Sec- 
ondly, by left-right symmetry {[Y^''\t)]^) = {[Yp^^\t)]'^) ; 
and they are obtained from the LE PU)) as 

{[Y^'-Ht)?) = {[Yp^Hm = 3NkBT/[4n^ki2p+inin) 

The third ingredient is that in terms of the mode am- 
plitudes Y^^\t) and Y^''\t) with the middle monomer 
fixed at TZ, the force (j)(t) on the middle monomer in Eq. 
(ID) is given by 



(/)(i) = 4k 



TT{2p+l) 

N 



Yp^'Ht)+Y(^\t)\. (12) 



With these ingredients, I consider an ensemble of phan- 
tom polymers (let us label it as ensemble X for future 
reference) in equilibrium at time 0~, such that for each 
of the polymers in this ensemble the velocity history of 
the middle monomer up to time t > tm is given by 



for some m and a sequence {Sri}, with to = 0. [One does 
not need not be alarmed that in Eq. ([T^ the velocity 
is represented by a sum over (5-f unctions in time: such 
a representation only simplifies the following analysis.] 
Note here, using Eq. (0), © and HI) that 



Yp{tJ) = Ypit^) - - 5n for ^ = 0, . . . , m. 



(14) 



Let us now determine the evolution of the forces the 
middle monomer of a phantom polymer, belonging to the 
ensemble X, experiences as a function of time. With the 
velocity history as in Eq. , using Eq. (jTU]) , the mode 
amplitudes for a member polymer of the ensemble evolve 
as 

Yp{t>t„,)^e-''-^'-'-^'^-Yp{tl) 



for Yp ^ [Y^\y^''\ Thereafter, using Eq. p^]fT5 



IpJt 



(j){t> tr. 



+4E 



{t-tm)/lp 



'^^'"+'^r^w,with 



N 



1 



+ - / rft'e-'=-(*-*')/^-[/('')(0 + 4'H0]- (16) 

Further, as I extend Eqs. (I14j|15p further back in time to 
write 

+ - dt'e-'=-(*'"-*')/>[/W(i') + /«(i')], (17) 



(t-t™)/7p 



Eq. (|T6| becomes 



N 



p 

.7r(2p+ 1) 



fcp(t-t„,_i)/7p 



p 

^^/g-fep(t-t')/7p[yW(i') + 4')(i')]- (18) 



IvJt 



Sri S{t - U) 



(^23) At this point, a comparison between Eqs. (|16l) and 

shows that the force (j){t > tm) on the middle monomer 



5 



of a polymer, belonging to ensemble X, can be extended 
all the way to time 0^, to write 



k 



E 



-kj,{t-ti)/-yp 



'4E^^^.i^°'W,with 



N 



+ i / dt'e-'=''(*-*')/T-[/W(i') + /i''(t')]-(19) 

It is now seen, by converting the sum to an integral, 
that 



8fc 



-kp{t~ti)/^p 



N 
2 5r, 



t-U 



-(t-U)/T 



(20) 



It is also seen, using Eq. PT|) . that {g{t))Q = 0, with the 
FDT 



N 



-(t-t')/r 



(21) 



assuming t > t' . Here (. . .)o denotes an average over the 
noise realizations, including an average over the configu- 
rations of the polymers of the ensemble X at t = 0^. 
The equilibrium condition for ensemble X at i = 0^ 
is necessary in Eqs. (|20II21I) to have {Yp^\t))Q- — 
{Yp^\t))Q- = by isotropy, and also to provide us with 

{[Y^''\mo- = {[Yp^'Hmo- =3NkBT/[4n^ki2p+iy] 
[these are the results of Eq. (ITU) ]. In other words, by 
means of Eqs. (|19ll2ip . I have derived the GLE ^ for 
the motion of the middle monomer of a phantom Rouse 
polymer, with 



(22) 



^{t)^2\l^ exp(-i/r). 



As claimed earlier, note that a comparison of Eq. ([22|) 
with Eq. (I^D|) shows that /i(t) is indeed built from the 
polymer's memory, i.e.. from the long relaxation times 
for fluctuations at long length-scales. In this context, 
it is worth recalling the work by Zwanzig [8], wherein 
the memory of the concerned particle arises due to its 
coupling to a thermal bath of harmonic oscillators, ef- 
fected via a Hamiltonian. For the present case, the mem- 
ory kernel ^,{t) arises due to the physical connectivity of 
the monomers: nevertheless, the fluctuation modes of the 



polymer effectively play the role of a thermal bath cou- 
pled to the motion of a monomer. 

Henceforth I drop the subscript '0' from the angular 
brackets; when applicable, an average over an equilibrium 
ensemble of polymers at t = will be understood. 



2. Anomalous dynamics for phantom polymers, and the 
derivation of Eq. for phantom polymers 

Having derived Eq. ([T]) for phantom polymers in Sec. 
Ill A 11 I now proceed to characterize the anomalous dy- 
namics for the middle monomer of a phantom polymer 
using Eq. Together, Eqs. ^ and (O yield 



jv{t) = dt' fi{t ~ t')v{t') + [g{t) + /^/2(t)] . (23) 

Jo V ' 



p(t) 



I then use Laplace transform to express 



V{t) = / dt'fi{t - t'Mt') + fN/2(t% (24) 

Jo 

where in the Laplace space /3{s)[j + jl{s)] — 1. Note that 
the noise term p{t) in Eq. ([23| also satisfies the FDT, 
as it should: the fN/2{t) part accounts for the viscous 
work involving 7 [this follows from the Rouse equation 
([3])[, while the g{t) accounts for fJ,{t), as in Eq. (|^^ . 

Before proceeding further I make the observation that 
Eqs. (|23ll24p are closed w.r.t. the velocity v{t) of the mid- 
dle monomer. In fact, Eqs. ()23II24|) provide the trajectory 
description in the (r, v) phase space, and clearly demon- 
strates that the stochastic process underlying the anoma- 
lous dynamics for phantom polymers is non-Markovian. 
I will return to this issue in Sec. IIVI 

The calculation for {v{t) ■ v{t')) based on Eq. ((M)) 
can be found in the Appendix. I quote the result below, 
assuming t > t': 



{v{t)-v{t')) =3kBTp{t~t'), 



(25) 



where /3{t) can be obtained by Laplace inverting the 
relation ^(s)[7 + p.{s)] = 1. Since from Eq. ([^ . 
/2(s) = 2y/TTjk{s + T^^)^^/'^, and therefore j3{s) = [7 -|- 
2\/Tryk(s-\-T~^)^^^^]^^ . Consequently, at long times, i.e., 
for times {t — t') ^ 7/fc but {t — t')<^T [this is possible 
in the limit of large N: recall that r = 77V^/(7r^fc)], the 
relation /3(s)[7 + /i(s)] = 1 in the Laplace space can be 
inverted to obtain 



{v{t)-v{t')) ^ (t-t')~'^'e-(*-*')/: 



(26) 



Subsequently, the result that the mean-square displace- 
ment (MSD) of the middle monomer increases ^ t^/"^ till 
time r and ^ t thereafter is obtained by integrating of 
{v{t) ■ v{t')) twice in time. 

It is important to note here that the following result 
transpires through the above exercise: if ^{t) ^ for 
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some a, then the MSD of the middle monomer will in- 
crease as t". 

Interestingly, instead of going through the Eqs. p3f 
I25p . one can choose an alternative route to arrive at Eq. 

that I took in Ref. I. Given the GLE ([1]), one arrives 
at Eq. ([2]), namely, at 

v{t) = - /" dt'a{t - t')${t') + h{t), (27) 

with jl{s)a{s) = 1 in the Laplace space, and the cor- 
responding FDT {h{t) ■ h{t'))o = {v{t) ■ w(t'))l^-=o = 
3kBT a{t — t'). Thereafter, in order to obtain the ve- 
locity autocorrelation function from Eq. (P7)l . one 
proceeds along the following manner. From the Rouse 
equation ([3]) or one needs to appreciate that 7/fc 
is proportional to the middle monomer's reaction time 



[For phantom polymers, the strain is due 

drn{t) 



to balance out the chain tensions 
drn{t) 



dn 



and 



dn 



n={N/2) + 

Thus, for (t — i') ^ 7/fc one is always 



n={N/2)- 

at a limit </> « 0, implying that 

{v{t) ■ v{t')) « {v{t) ■ v{t'))$=o = iksTait - t'), (28) 

i.e., I arrive back at Eq. ([^5]) [9], which I have elaborated 
via Eqs. pll25l) . 

To summarize Sec. |ll]so far: (i) up to Eq. ([26|) I have 
provided an elaborate derivation for the GLE formula- 
tion for the anomalous dynamics of phantom polymers, 
explicitly demonstrating that the stochastic process un- 
derlying the anomalous dynamics of the middle monomer 
is non-Markovian. (ii) Given that the GLE formulation 
([T][2]) dictates that if /i(i) ^ then the anomalous 

dynamics exponent is also a, characterization of fi(t) is 
sufficient to lead one to the anomalous dynamical behav- 
ior of phantom polymers. 

In the following subsections III BMII Dl and in Sec. 
mil I will consider the anomalous dynamics for self- 
avoiding Rouse and Zimm polymers, polymer transloca- 
tion through a narrow pore in a membrane, polymer rep- 
tation and polymer melts. These systems are not analyt- 
ically tractable like phantom polymers have been. Never- 
theless, as I have demonstrated above for phantom poly- 
mers, and summarized in Ref. I, characterizing /i(i) will 
turn out to be powerful enough to arrive at the anoma- 
lous dynamical behavior for all these polymeric systems 
via the GLE description ([TJ2]). To this end, one needs to 
identify the generic properties involving fj,(t) and FDT; 
these are carried out in Sec. Ill A 31 below. 



3. Three generic issues for polymeric systems: (a) fi{t) as 
the polymer's mean relaxation response to local strain, (b) 
fluctuation- dissipation theorem, and (c) Eq. \23\) for 
polymeric systems in general 

I begin by having observed the following from Eq. 
([T]): nit) is the mean response of the polymer to a lo- 



cal strain 

to the alteration in the chain tensions k 



dn 



and k 



dn 



--(N/2) + 

caused by moving the middle 



n={N/2)- 

monomer by a distance Sr at time to and fixing it at its 
new position Vt. Then on average, the local strain then 

relaxes in time ~ (t — io)^"'"^^, i-e., (0(0 ■ 4'{to))v=o ~ 
{t - to)-^/^.] Based on this observation, I now discuss 
three generic issues that would be instrumental for arriv- 
ing at the anomalous dynamics for polymeric systems, as 
these systems are not analytically tractable like phantom 
polymers have been. 

(a) While the mean response of a polymer to a lo- 
cal strain depends on how the strain is created, 
the identification of ii{t) as the polymer's mean lo- 
cal strain relaxation response alone allows one to 
write down the GLE, as I argue below. First, 
given the identification of fj,{t) as the polymer's 
mean local strain relaxation response one can al- 
ways write the stochastic Eq. ([l} with {g{t)) = 0, 
which holds by definition. Next, to obtain the 
FDT, consider Eq. ([1]) for an ensemble of poly- 
mers with v(t) = Vt and (t>{to) = go, a specific 
value. For such an ensemble g{t) = (t>{t), and since 
IJi{t) is the polymer's mean local strain relaxation 
response, {4>{t) ■ cf){to)) — ffo /^(t — to) for t > to- 
Extending this to the dynamics of a polymer in an 
equilibrium ensemble (where go is also chosen from 
the equilibrium ensemble), one has {g{t) ■ g{to)) = 

{${t) ■ ${to))y=o = {^^{t))v=0 f^{t - to). 

(b) Conversely — i.e., turning the argument of (a) 
around — the existence of the FDT {g{t) ■ g{to)) = 
C ii{t — to) for some constant C, along with {g{t)) — 
also implies that the polymer's mean relaxation 
response to local strains — created at time to — 
is cx /i(t — to). Let us assume that the poly- 
mer's mean relaxation response to local strains, 
created at time to, is given hy rj{t — to) for some 
function -q. Using (a) above, one deduces that 
the corresponding FDT will have to take the form 
(5(0 -5(^0)) = {g^{to))r]{t~'to), which, by construc- 
tion equals Cfj,{t — to); i.e., r]{t) oc fJ-(t). 

(c) In addition, for a tagged monomer in a polymeric 
system, one anticipates Eq. ([SJ to generically hold 
true, with 7 being an effective damping coefficient 
of the surrounding medium (that includes the in- 
teractions with other monomers!). Together, they 
then yield Eq. which is a closed form equa- 
tion for v{t). It is Eq. (HJ) that tells us that if 
the monomer makes a move at any time, there is 
an enhanced chance to undo this move in subse- 
quent times; this is where anomalous dynamics in 
polymeric systems stem from. 
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In Sees. HTBlllTDj andlllll in order to arrive at the GLE 
formulation (HJ^]) for the anomalous dynamics for self- 
avoiding Rouse and Zimm polymers, polymer translo- 
cation through a narrow pore in a membrane, polymer 
reptation and polymer melts, I will use (a-c) above, since 
these systems are not analytically tractable like phan- 
tom polymers have been. In some cases I will obtain fi(t) 
directly, and in others I will calculate {(f>{t) ■ (j>{t'))^^Q. 
Wherever possible, I will do both. 

B. The GLE formulation ([l][2]) for the anomalous 
dynamics of self-avoiding Rouse and Zimm polymers 

The monomers of a self-avoiding polymer interact over 
a long-range, which prohibits one from writing down an 
exact equation for the velocities of the monomers in terms 
of the forces they experience. However, some properties 
of self-avoiding polymers are well-known; two of them I 
need here for a polymer of length N are: (i) the ter- 
minal time T scales ~ N^+'^'^ for a Rouse polymer, and 
as ~ N^'^ for a Zimm polymer and (ii) the entropic 
spring constant of a polymer scales as N~^'^ Here 
1/ is the Flory exponent, in 3D i/ « 0.588, and in 2D 
V = 3/4. With these properties I now characterize fi(t) 
as the polymers' mean local strain relaxation response, 
as described in Sec. Ill A 31 
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FIG. 1: (color online) Simulation data for fi{t) for a self- 
avoiding Rouse polymer in a Monte Carlo lattice polymer 
model for iV = 400, obtained from Eq. ^ and the FDT 
(see text for details). With 2i^/(l + 2!/) ^ 0.54 in 3D the solid 
(black) curve corresponds to fi{t) = 0.275f~'''^^ exp(— f/r), 
where r was measured to be ~ 223801. Figure reproduced 
from Ref. [lol | with permission from American Institute of 
Physics. 

Imagine that one moves the middle monomer of a self- 
avoiding polymer by a small distance 6r at t — and 



holds it at its new position \/t > 0, corresponding to 
v{t) — 5rS{t). Following the terminal time scaling for 
self-avoiding polymers (i), at time t, counting away from 
the middle monomer, all the monomers within a back- 
bone distance nt ~ t^/(^+'^'^) for a Rouse, and ~ t^/(^'^) 
for a Zimm polymer equilibrate to the new position of 
the middle monomer. However, since the rest (N — nt) 
monomers are not equilibrated to the new position of 
the middle monomer at time t, these nt monomers are 
stretched by a distance 6r. With the entropic spring 
constant of these nt equilibrated monomers scaling ^ 
n^^" [following property (ii)], the mean force the mid- 
dle monomer will experience at its new position is given 
by ${t) - nt^'^i-S'r) - t-2^/(i+2'^)(-fr) for a Rouse, 
and 0(<) ^ n^"^^ {—5r) ^ t^'^/'^{~5r) for a Zimm polymer 
[force ~ (spring constant) x (stretching distance)[. This 
power-law behavior lasts only till the terminal time r. 
[for a phantom polymer, the time behavior of Eq. is 
recovered from this of argument upon simply replacing v 
by 1/2.] In the light of point (a) in Sec. Ill A 31 the above 
leads one to the GLE formulation ([T]) for self-avoiding 
Rouse and Zimm polymers, with the corresponding FDT. 



Further, having anticipated, as in point (c) in Sec. 
Ill A 31 that an effective Eq. does hold true, I demon- 
strate the FDT directly for self-avoiding Rouse polymers 
in the following manner. Imagine that I hold fixed the 
middle monomer of a self-avoiding Rouse polymer, and 
keep taking snapshots of the polymer configuration at 
fixed intervals of time, i.e., at t = io, (io + At), (to + 
2 At),.... Afterwards, I take the snapshot at time t 
and evolve it K number of times over a short time in- 
terval 5t and note down the average displacement 5x 
over these K evolutions. A look at Eq. (pS)) then tells 
us that this averaging process kills the term, and 

since v{t) = by construction, the result is a quan- 
tity 5x, which, in the limit of [K — >■ oo, 6t 0), is 

simply proportional to the average of g{t). These Sx 
values for t — to, (to + At), (to + 2 At), . . . can then be 
used to calculate /i(t — t') = {g{t) ■ g{t'))ij=o, proxied by 

{dx(t) ■ 5x{t')). Such an exercise was indeed performed 
m Ref. [ICi]: using a Monte Carlo lattice polymer model 
with N = 400, At = 50, if = 10^ and 5t ^ I (see 
Ref. for model details). The corresponding result, 
confirming the ^-21^/(1+21/) power-law behavior of /i(t) is 
shown in Fig. [T] 



Given that if /i(t) ^ t " then the anomalous dynamics 
exponent is also a, the well-known scaling of the MSD 
~ i(2^')/(i-H2i') foj. self-avoiding Rouse and ~ t^/^ for self- 
avoiding Zimm polymers up to time t and ^ t thereafter 
is trivially reproduced from the GLE formulation ([TJ2]) 

©[HI . 
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FIG. 2: (color online) The mean relaxation response ~ local strain created by injecting n monomers at 

t = into a self- avoiding Rouse polymer tethered on a membrane, at the tether point. The data are obtained using a Monte 
Carlo lattice polymer model, (a) in 2D: for — 250, n = 5. In 2D (1 + + 2u) — 0.7; the solid line corresponds to a 

power-law (b) in 3D: for n = 10 and iV = 50 (red), 100 (green) and 150 (blue). In 3D (1 -I- i/)/(l -I- 2v) « 0.73; the solid 

line corresponds to a power-law t""'^^. The choice of the scaling variable also confirms the terminal exponential decay of ij,{t) 
behaving ~ exp(— t/r), with r ~ N^^'^" for self-avoiding Rouse polymers. The figures are taken from Refs. and [l5j, where 
the details of the polymer model and simulation details can be found, with permissions from Institute of Physics Publishing 
Ltd., UK (original articles published on 25 January 2008 and 8 October 2007 respectively). 



C. The GLE formulation (IT][2]) for the anomalous 
dynamics of unbiased polymer translocation through 
a narrow pore in a membrane 

Polymer translocation is a process where a polymer 
passes through a narrow pore in a membrane. Of interest 
here is the so-called unbiased (i.e., in the absence of any 
force or field) translocation: the polymer passes through 
the pore purely due to thermal fluctuations, and the dy- 
namics is anomalous For this system, below I now 
characterize fi(t) as the polymers' mean local strain re- 
laxation response, as described in point (a) of Sec. Ill A 31 
in the following manner. 

A translocating polymer consists of two sub-polymers 
— one on each side of the membrane — exchanging 
monomers through the pore. When a monomer translo- 
cates, the polymer locally stretches on the side the 
monomer translocates from, and locally compresses on 
the other. Consequently, the polymer's chain tensions 
at the pore changes: it increases on the side of the 
membrane which the monomer translocates from, and 
decreases on the other. The relevant mean polymeric 
response therefore, is to a (local) strain due to extra 
monomers injection at the tether point of a polymer teth- 
ered on a membrane. 

Consider the case where n extra monomers are injected 
at t = into a polymer tethered on a membrane, at 
the tether point. First, it has been shown in Ref. (T3| 
that for phantom Rouse polymers the mean response to 
such a strain is given by ^(i) ~ t~^e~*^'^ , with r ~ N^. 



For self-avoiding polymers /i(i) is obtained as follows. 
Given that the terminal time r for a tethered polymer 
of length N scale N^+^'' for a Rouse and N^" 
for a Zimm polymer, at time t, counting away from the 
tether point, all the monomers within a backbone dis- 
tance nt ^ ti/(i+2'') for a Rouse, and t'^/i^'^) for a 
Zimm polymer, equilibrate to the injected monomers. 
The real space extent of rit monomers is r(nj) ~ n^, 
but since the rest {N — rit) monomers are not equili- 
brated to the injected monomers at time t, there are 
(rit + n) monomers squeezed in a space that extends 
only to r{nt). The corresponding compressive force [force 
— (spring constant) x (stretching distance)] from these 
(nt -|- n) monomers, felt at the pore, and hence fi(t), 
is the given by ~ n^'^'^[Sr{nt)] ^ n^'^'^ n[dr{nt) / dnt] — 
vnn'^^^^'^\ which scales ^ t-0-+^)/('^+'^^^) for a Rouse and 
~ ^-(i+i'j/CSi') fQj. Zimm polymer. (Once again, this 
behavior lasts only till the terminal time r.) The mean 
relaxation response to the local strain for a self-avoiding 
Rouse polymer in 2 and 3D using a Monte Carlo lattice 
polymer model — originally explicitly evaluated in Refs. 
[11 [ill — are shown in Fig. [H (see Refs. [H ^ for the 
model details). 

The above thus implies that the anomalous dynam- 
ics for polymer translocation is described by the GLE, 
resuhing in the scaling of the MSD - t(i+'')/(i+2'') for 
self-avoiding Rouse and ^ tO-+'^)/i^^) for self-avoiding 
Zimm polymers up to time r and ^ t thereafter. Conse- 
quently, the pore-blockade time scales ^ N"^ for aphan- 
tom Rouse ~ N'^+'' for self-avoiding Rouse (l4.fT5t. 
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and iV^+^'^ for self-avoiding Zimm (Tsl . [T3| polymers. 

D. The GLE formulation ([l][2]) for a reptating 
polymer 

The notion of polymer reptation was originally intro- 
duced by De Gennes in the context of gel electrophore- 
sis Gel electrophoresis is a method to isolate DNA 
fragments. In a typical experiment, agarose powder is 
dissolved in water, and upon waiting for a sufficient time 
the agarose forms the cross- links of a gel. The sizes of the 
pores within the gel is controlled by the concentration of 
agarose. When DNA molecules are injected into this gel, 
they collect their charge from the solution, and with the 
application of an electric field — typically of strengths 
of a few volts per cm — the DNA molecules are driven 
through the pores in the gel. For a given pore size of 
the gel, the mobility of the DNA molecules caused by 
the electric field crucially depends on their lengths, lead- 
ing to length-dependent segregation, which in turn allows 
one to separate the DNA strands by their lengths. The 
DNA persistent length under a typical gel electrophoresis 
experimental condition is in the range 400-800 A, which 
is typically the size of the pores in the gel itself. In other 
words, the gel prevents the polymer to move transverse 
to its configuration, and only the longitudinal, or curvi- 
linear, motion of the DNA contributes to its motion un- 
der the applied electric field — this is the central idea of 
polymer reptation. 




FIG. 3: A two-dimensional lattice representation of the repton 
model. The pores are represented by cells formed by a square 
lattice, while the polymer is represented by points ("reptons") 
connected by bonds. The rules governing this model are de- 
scribed in the text. 

I consider here the three-dimensional repton model, 
one of the simplest models with explicit reptation (lon- 



gitudinal/curvilinear) moves. The model was originally 
conceived by Rubinstein [l^ to describe the motion of 
entangled polymers, and was later co-opted by Duke [l^l 
to describe the mobility of DNA polymers in gel elec- 
trophoresis. A two-dimensional lattice representation of 
this model is shown in Fig. |31 the network of the agarose 
gel is represented by the crossing points of the grid lines of 
the lattice, whereas the cells of the lattice correspond to 
the gel pores. The DNA polymer, represented by points 
connected by bonds, resides in the gel. The points are 
called "reptons": along the backbone of the polymer two 
consecutive reptons are one persistent length apart. This 
allows more than one repton to occupy the same cell, but 
simultaneously two neighboring reptons cannot be more 
than one cell apart. The movements of the reptons are 
given by the following dynamical rules. A repton in the 
interior of the polymer can move to one of its adjacent 
cells, provided one of its neighbor reptons is already in 
that cell, and the other neighbor repton is in the cell it 
leaves. The two end reptons can move into any of the 
neighboring cells, so long as the bond between the end 
repton and its neighbor is separated by not more than 
one cell. 

The dynamical rules of the repton model ensures that 
while the transverse motion for a repton in the interior 
of the polymer remains blocked, it can only move to- 
wards (and away from) its neighbors: this encodes ex- 
plicitly the longitudinal or curvilinear motion of the poly- 
mer in the repton model. Further, in any spatial dimen- 
sion, the curvilinear conformation of a reptating polymer 
can be described by a string of Os and Is, correspond- 
ing to whether two consecutive monomers reside within 
the same cell or not — this forms the "curvilinear co- 
ordinate" description of the polymer. (E.g., the curvilin- 
ear co-ordinate description of the two-dimensional poly- 
mer in Fig. [21 from the left end to the right, is given 
by the sequence 01111101101111110111111011101). The 
Os in the curvilinear co-ordinate description of the repton 
model is also known as the "stored lengths", as they corre- 
spond to local compression of the polymer: the number of 
stored lengths within any given lattice site is clearly one 
less than the number of monomers occupying that site. 
The dynamics of the repton model is then summarized by 
the two following simple rules: within the polymer mo- 
tion takes place only by exchanging a with a 1, while 
at the ends of the polymer, a zero can become a 1, or 
vice versa. The equilibrium density of stored lengths is a 
parameter in the repton model, for the simulation results 
reported in this paper, this parameter was chosen to be 
1/3. 

The dynamics of reptation is anomalous 0; i.e., for 
a polymer of N reptons the MSD of the middle repton 
increases ~ t^/'* up to time ~ N^; thereafter, the MSD 
for the middle monomer increases ^ i^/^ till time ~ N^, 
after which the motion becomes difi^usive. I confirm this 
anomalous dynamical behavior of the middle repton in 
Fig. m By the time the MSD behavior becomes diffu- 
sive, the middle monomer displaces itself by the size of 
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the polymer, and since by construction the repton model 
describes a phantom polymer (i.e., in equilibrium the size 
of a reptating polymer of N reptons scales ~ JV^/^, the 
diffusion coefficient of the polymer scales as (2ll . [2^ . 
In fact, given that the polymer's motion is only curvilin- 
ear, and the polymer's contour itself performs a random 
walk in space, one realizes that if this anomalous dynam- 
ics of the middle monomer is translated back into curvi- 
linear co-ordinates, then the MSD of the middle repton 
increases ^ t^/^ up to time ~ A^^, and becomes diffusive 
thereafter. It is this anomalous dynamics of the middle 
monomer in curvilinear co-ordinates that I show to stem 
from the GLE ([T][2]). In particular, I demonstrate that 
the forces on the middle monomer in the curvilinear co- 
ordinate, at zero middle repton velocity, satisfy the FDT, 
leading one to fi{t) ^ exp(— i/r) with r ~ iV^. As 
described in point (b) of Sec. Ill A 31 this establishes the 
GLE formulation for the repton model. 
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FIG. 4: The MSD (Ar^(t)) of the middle repton for polymers 
with iV reptons: TV = 201 (circles), 501 (pluses) and 1001 
(crosses). Data averaged over 1,024 independent polymer re- 
alizations. Upto time 7=^ 0.25iV^ the MSD behaves as t^^^, 
after which it crosses over to ~ i}^"^ behavior. The onset of 
the diffusive behavior at time ~ A'^'^ is not shown, although 
increased slopes at late times for shorter polymers can be seen 
clearly. 

The above implies that the anomalous dynamics for 
polymer reptation is described by the GLE, resulting in 
the scaling of the MSD in the curvilinear co-ordinate 
~ t^l"^ till time r and diffusive thereafter, as confirmed 
in Fig. 13) In the curvilinear co-ordinate, the polymeric 
motion is described by the motion of the stored lengths 
or the Is, which undergo the so-called tagged-particle dif- 
fusion |22l,[2l|. At zero velocity of the middle repton, the 
force it experiences is proportional to the curvilinear gra- 
dient of stored lengths at its location. In order to obtain 
the correlation in the forces on the middle repton at zero 
velocity, I number the reptons 1, 2, . . . , iV, hold the mid- 
dle repton (tagged by fixed, and monitor the differ- 
ence g{t) = PN/2,5(t)- PN/2-5{t), where Pm,n{t) is the of 



number of stored lengths between reptons m and m + n. 
The correlation in the force on the middle repton at zero 
velocity is then proxied by {g{t)g{t')). The expected be- 
havior of p{t - t') oc {g{t)g{t')) - i-i/2gxp(_i/r) is 
shown in Fig. [SJ with r ~ 7V^. 
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FIG. 5: The ~ f"^/^ exp(-f/r) behavior of ii{t), obtained 
via the FDT in the curvilinear co-ordinate for the anomalous 
dynamics of the middle repton for three different values of 
N: N = 201 (circles), 501 (pluses) and 1001 (crosses). Data 
averaged over 1,024 independent polymer realizations. The 
solid line corresponds to the power-law t~^^'^ . The choice of 
the scaling variable also confirms the terminal exponential 
decay of p{t) behaving ~ exp(— f/r), with t ~ A''^. 



III. THE GLE FORMULATION ([T][2]) FOR THE 
ANOMALOUS DYNAMICS OF POLYMER 
MELTS 

A. Polymer melts and reptation theory 

Mobility of individual polymers in systems containing 
long polymers decreases with increasing density, caus- 
ing their dynamics to slow down. A notable manifesta- 
tion of slow dynamics in dense polymeric systems such 
as polymer melts is that their viscosity scales with their 
molecular weight as a power law, with the viscosity ex- 
ponent m^; experimentally measured to be 3.4 ±0.2 over 
an impressive range of molecular weights and chemical 
compositions [23|. The fact that the slow dynamics of 
dense polymeric systems can actually be attributed to 
the polymers becoming entangled, especially if tliepoly- 
mers are long, was the seminal idea of De Gennes [1|. He 
argued that entanglement restricts the polymers to only 
sliding past each other — thereby reducing their mobility 
— as it is impossible for them to slide across each other. 
Indeed, the dynamics of a melt can be quantitatively un- 
derstood by taking the entanglement idea to analyze the 
motion of a single polymer in a melt P, Q — referred to 
as the tagged polymer henceforth. 
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To approach the scahng of by obtaining that of D, 
De Gennes considered the meU to only represent a rigid 
static network for the tagged polymer, with the tagged 
polymer reptating through this network. Such a formu- 
lation reduces the many-polymer problem of a melt to 
an effective single (tagged) polymer problem, wherein the 
movements of the tagged polymer is restricted only to 
reptation, i.e., stored length transportation due to (lon- 
gitudinal) fluctuations along its contour. By using the 
longest available time-scale ~ iV^ for the stored length 
distribution within the tagged polymer's contour in this 
reduced problem, he obtained D ~ N~'^, i.e., to,, = 3 
P, 0, [l^ . [The discrepancy between the predictions of 
reptation theory and the experimental value for m,, re- 
mains to this day, and is beyond the scope of this paper.] 




FIG. 6: (color online) Illustration of the two-dimensional ver- 
sion of the lattice polymer model. Polymers are shown by 
darker colors and their contours in lighter colors. In the up- 
per polymer, interior monomers 2, 4, 6, 9, 10 and 11 can either 
move along the contour, or move sideways; monomer 7 can 
join either 6 or 8; the end monomers 1 and 12 can move to 
any empty nearest-neighbor site. In the lower polymer, inte- 
rior monomers 3, 5, 6, 10 and 11 can either move along the 
contour, or move sideways; monomer 1 can move to any empty 
nearest-neighbor site, and monomer 12 can join its neighbor 
11. In this configuration, because of the self- and mutually- 
avoiding property of the polymers, all other monomers can- 
not make a move. Statistically once per unit of time, each 
monomer attempts to move up and down along the contour, 
as well as sideways. 

The prediction of reptation theory for the anomalous 
dynamics of the middle monomer of the tagged polymer 
in a melt, beyond th so-called entanglement time Tg is 
the same as that for the repton model discussed in Sec. 
IIIDI at times ^ t^, the MSD of the middle monomer 
increases ^ t^/** up to time ^ iV^; thereafter, it increases 
~ t^^"^ till time ^ N^, after which the motion becomes 
diffusive. Such a picture is consistent with the fact that 
by the time the dynamics becomes diffusive, the mid- 



dle monomer displaces itself by the size of the polymer, 
which scales as ^ N^^^ in a melt [i], yielding the result 
D ~ iV^^. For the GLE description of the anomalous 
dynamics for a many-body problem like a polymer melt, 
I make use of computer simulations to demonstrate the 
FDT for the force experienced by the middle monomer 
of the tagged polymer, as described below. 



B. Simulation details 

The simulations have been performed with the lattice 
polymer model described in Ref. [25j. This model com- 
bines a very high computational efficiency with realis- 
tic polymer dynamics. All the simulations are three- 
dimensional, with the lattice polymers residing on a face- 
centered-cubic lattice. Polymer contours are self- and 
mutually-avoiding. Monomers adjacent in the string are 
located either in the same, or in neighboring lattice sites. 
Multiple occupation of lattice sites is not allowed, ex- 
cept for a string of adjacent monomers belonging to the 
same polymer. A two-dimensional version of the model 
is illustrated in Fig. [S) 




FIG. 7: fi{t) for the middle monomer of a tagged polymer in 
a polymer melt: = 750 (circles), 1500 (pluses). Data av- 
eraged over 512 independent polymer melt realizations. The 
solid line corresponds to the power-law t~^^*. Beyond the 
entanglement time Te, the data are consistent with an initial 
power-law behavior fj,{t) ~ t~ 
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The simulations are performed for a system of size 60'^ 
with an overall monomer density unity per lattice site. 
The polymers move through a sequence of random single- 
monomer hops to neighboring lattice sites. These hops 
can be along the contour of the polymer, thus explic- 
itly providing reptation dynamics. They can also change 
the contour "sideways", providing Rouse dynamics (Fig. 
[6]). Each monomer attempts to move along the contour, 
as well as sideways stretching or reducing the backbone, 
statistically once per unit of time. Due to the possibility 
that adjacent monomers belonging to the same polymer 
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polymeric system 






self-avoiding Rouse 


~ t-2-/(i+2-) exp(-t/r) 


~ t-i/(i+2>^)^ tiU T and ~ F thereafter 
{[x{t) ~ x{0)]) ~ t'^'^/iW,^)F till r, and ~ tF thereafter 


self-avoiding Zinini 


~ t"^''^exp(-t/r) 


~ t'^/^F tiU r and ~ F thereafter 
{[x{i) - f(0)]) ~ t^^^F till T and ~ tF thereafter (not shown here) 


reptation; repton model 
(curvilinear co-ordinate) 


~ t'^^' exp{-t/T) 


~ t^^^'^ till r and ~ t thereafter 
{[x{t) - x{Q)]) ~ t^/^F till r and t ~ F thereafter 


melt (reptation theory) 


~ between Te and ~ A'^^ 
(as much as data can resolve) 


~ t'^'^ between Te and ~ iV^ (*) 
{[x(t) — x{Q)]) ~ t^^'^F between and ~ A''^ (not shown here) 



polymeric system 




number of translocated monomers in time t 


phantom Rouse 
field-driven polymer translocation (3D) 


~t"^exp(-t/r) 
(r not translocation time) 


~ t 

translocation time ~ A'^ (not shown here) 


self-avoiding Rouse 
field-driven polymer translocation (3D) 


^ i-(i+'^)/(i+2.) exp(-Vr) 
(r not translocation time) 


^ j(l+i/)/{l-l-2i.) 

translocation time Af(^+2'"'''(^+''' 


self-avoiding Zimm 
field-driven polymer translocation 


_ ^-(i+.)/{3.) exp(-t/r) 
(r not translocation time) 


translocation time ~ j\r^''/(i+'') (not shown here) 


self-avoiding Rouse 
pulled polymer translocation 


~ t~^/^exp(-f/r) 
r ~ A*'^ (r not translocation time) 


,1/2 

translocation time ~ A''^ 



polymeric system 


flit) 


number of monomers adsorbed in time t 


phantom Rouse 
polymer adsorption (3D) 


~ exp(— f/r) 
(r not adsorption time) 


~ t 

adsorption time ~ N (not shown here) 


self-avoiding Rouse 
polymer adsorption (3D) 


_^-(i+.)/(i+2,.) exp(-t/r) 
(r not adsorption time) 


^ i(l+'^)/{l+2i^) 

adsorption time ~ Ar(i+2'')/(i+'') 


self-avoiding Zimm 
polymer adsorption 


_ j-(l + .)/{3.)g^p(_^/^-) 

(r not adsorption time) 


^ ^{l+"-)/(3--) 

adsorption time ~ (not shown here) 



TABLE II: The expected drift results, from the analysis in Sees. IIV Al and IIVBI — when the systems start evolving at t = 
under the application of weak forces — following the results of Table [T] The top table relates to Sec. IIV Al and the bottom 
two tables relate to Sec. IIV Bl Note in Table [III that if /i(f) ^ t~" , then the drift exponent is a. Note also that not all 
results are explicitly demonstrated in Sees. IIV Al and IIV Bl using computer simulations; these are clearly marked. (*) Te is the 
entanglement time for polymer melts within the reptation theory Q. 



can occupy the same site, overall approximately 40% of 
the sites typically remain empty. The number of stored 
lengths within any given lattice site is one less than the 
number of monomers occupying that site. 

Initial thermalizations were performed as follows: com- 
pletely crumpled up polymers are placed in lattice sites 
at random. The system is then brought to equilibrium 
by letting it evolve for Teq(iV) units of time, with a com- 
bination of intermediate redistribution of stored lengths 
within a polymer [2^ ; for the polymer lengths concerned 
hereToq(A^) > 10^. The entanglement time for this model 
is w 10^. 



C. The FDT for the forces on the middle monomer 
of the tagged polymer in a polymer melt 

First, I anticipate, as in point (c) in Sec. Ill A 3[ that an 
effective Eq. ((231) does hold true for the middle monomer 



of the tagged polymer in a polymer melt. Then I follow 
the same procedure along the as I did for demonstrating 
the FDT for the case of self-avoiding Rouse polymers; 
namely, I hold the middle monomer fixed and keep tak- 
ing snapshots of the polymer configuration at fixed in- 
tervals of time, i.e., at t = to, (^o + At), (to + 2 At), . . .. 
Afterwards, I take the snapshot at time t and evolve it 
K number of times over a short time interval St and note 
down the average displacement Sx over these K evolu- 
tions. Once again, a look at Eq. tells us that this 
averaging process kills the fN/2 term, and since v{t) = 

by construction, the result is a quantity Sx, which, in the 
limit of {K oo, (5t 0), is proportional to the average 

of g{t). These S^x values for t = to, (to+At), (to + 2At), . . . 
can then be used to calculate /i(t — t') = {g{t) ■ g{t'))fj=o, 

proxied by {Sx{t) ■ Sx{t')). 

The corresponding result for N = 750 and 1500 for 
At = 10, 000 and 50, 000 respectively, with K = 100, 000 
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and (5i = 1 is shown in Fig. [T] Beyond the entanglement 
time Tg , the data are consistent with an initial power-law 
behavior ^ (but are not sufficient to conclude 

what follows afterwards). Nevertheless, this exercise pro- 
vides a clear indication that the anomalous dynamics of 
the middle monomer of the tagged polymer is described 
by the GLE between the entanglement time Tg and ^ ^ . 



IV. DRIFTS IN POLYMERIC SYSTEMS 
DRIVEN BY WEAK FORCES: ANALOGUE OF 
THE NERNST-EINSTEIN RELATION 

Until now, I have considered the classical polymeric 
systems with anomalous dynamics in the absence of 
external driving forces, and demonstrated that their 
anomalous dynamics are described by a unified GLE 
scheme (IT][2|) with power-law memory kernels. I now 
demonstrate that their GLE formulation is robust: it also 
adequately describes drifts in polymeric systems driven 
by weak forces, as it should. 
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FIG. 8: Confirmation of Eq. (|32p for an end monomer of a 
self-avoiding Rouse polymer: data collapse for the displace- 
ment {Az(f)) in the f-direction, under the application of a 
force F — Fz from t = onwards. Data are obtained using 
the same Monte Carlo lattice polymer model that was used 
to generate Fig. [T] iV = 300, FA/(fcsT) = 0.075 (circles) and 
N = iOO, FX/iksT) = 0.1 (pluses), where A is the lattice 
constant. Data averaged over 8,192 equilibrated polymer re- 
alizations at t = 0. Following Eq. p2p . {Az{t)) is expected 
to behave as t2f/{i+2i^)^ ^ t°''^^F till the terminal relaxation 
time r ^ N^^^" , after which it should behave as tF. Note 
the labeling of the axes: they confirm that in the steady state 
{t > r) the drift scales as 1/A'^ as it should. 

The expected drift results — when the polymeric sys- 
tem starts to evolve at t = under the application of 
external forces — are shown in Table |lTl Note in Table 
|TT]that if iJi{t) ~ then the drift exponent is also a: 1 
elaborate on this below. Note also that not all results in 



Table HIl are demonstrated in Sees. IIV Al and IIV Bl using 
computer simulations; these are clearly marked. 



A. The standard extension of the GLE formulation 
to driven polymeric systems 

The external constant force F is switched on a specified 
monomer over an equilibrated polymer ensemble a.t t — 
0, and is assumed to be weak enough that it does not 
distort the memory kernel [or the FDT for g{t)\. With 
these assumptions — if the system's dynamics is not rate- 
limited by some other process — Eqs. (IT][5]) retain their 
forms, however, Eq. needs to be modified to include 
the external force, so that it becomes 



(29) 



from where I have dropped the subscript N/2 from f{t). 
Together with Eqs. (PH]), Eq. ^ provides us with 



v{t) 



dt' p(t-t')[g{t') + F + f{t% 



(30) 



where /3(i) is defined in Eq. (P^ . 
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FIG. 9: Confirmation of Eq. (|32p for the middle repton of the 
repton model. Data collapse for the curvilinear displacement 
(As(t)), under the application of a curvilinear force F from 
t = onwards: TV = 201,FA/(fcsr) = 0.1 (circles), TV = 
501,FA/(fesT) = 0.25 (pluses) and TV = 1001, FA/(fcsr) = 
0.5 (crosses), where A is the lattice constant. Data averaged 
over 16,384 independent polymer realizations at equilibrium 
at f = 0. Following Eq. ([32} {^s{t)) is expected to behave 
as t^^^F till the terminal relaxation time r ~ TV^, after which 
it should behave as tF. Note the labeling of the axes: they 
confirm that in the steady state {t > r) the drift scales as 
1/TV as it should. 

An ensemble average of Eq. (pO)) implies that 



(vit)) = / dt' /3{t-t')F. 
Jo 



(31) 
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With ^ t "exp(— t/r) for some a, at long times 
(times still < t), Eq. ([211) yields 

(t7(i)) ^ r-iF; i.e., {[x{t) ^ x{0)]) - , (32) 

which lasts till time r, after which {x{t)) increases as tF. 
Such behavior can be thought of as the analogue of the 
Nernst-Einstein relation [6]. 

I now demonstrate Eq. for a self-avoiding Rouse 
polymer, and the repton model in the curvilinear co- 
ordinate. 



1. A self- avoiding Rouse polymer under a weak force on an 
end monomer 

Using the Monte Carlo lattice polymer model that was 
used to generate Fig. [1] (details of the model can be 
found in Ref. [13 )j a force F = Fz was applied at t = 
onwards on one of the end monomers of a set of equili- 
brated self-avoiding Rouse polymer realizations ai t — Q, 
and the average displacement (Az(t)) of this monomer 
was tracked as a function of time. The result, presented 
in Fig. [51 confirms Eq. ((5^ . Note the labeling of the 
axes in Fig. (5) they confirm that in the steady state 
{t > r) the drift scales as 1/A^ as it should. 

2. The repton model with a weak curvilinear force on the 
middle repton 

Similarly, for the repton model a curvilinear force F 
was applied from t = onwards on the middle repton 
of a set of equilibrated polymer realizations to t = 0, 
and the average curvilinear displacement (As(t)) of this 
monomer was tracked as a function of time. The result, 
presented in Fig. [^1 again confirms Eq. (P^ . Note the 
labeling of the axes in Fig. [HI they confirm that in the 
steady state {t > r) the drift scales as 1/A^ as it should. 

B. A different extension of the GLE description to 
driven polymeric systems 

If the external force is not acting on a specified 
monomer (or the system's dynamics is rate-limited by 
some other process), then a different generalization of the 
GLE ([T][5]) is required to obtain the drifts. This holds for 
dynamics of field-driven translocation in 3D f2f\, poly- 
mer adsorption on a solid surface in 3D [28\ , and poly- 
mer translocation by means of a pulling force F (30| . For 
the first two cases, the force does not act on a specific 
monomer; instead, it acts respectively on the monomer 
that is in the pore, and on the monomer that is just in 
contact with the adsorbing surface. However, for poly- 
mer translocation by means of a pulling force, the force 
does act on an end monomer, but the translocation dy- 
namics is rate-limited by polymer dynamics at the pore. 



For these problems, a description based on Eq. ^ alone 
is sufficient, with (/)(0) ^ 0, provided, once again, that 
the force F is small enough that it does not distort the 
memory kernel or the FDT for g{t). In that case, Eq. ([T]) 
reads 

0(t) = 0(0) - / dt'^l{t - t')v{t') + (33) 

Jo 

leading one to 

m) = f dt' a{t - 1') am - ${t')]), (34) 

Jo 

which is essentially an extension of Eq. If {(f>{t)) be- 
comes a constant fairly soon after the force starts acting 
on the polymer, then from Eq. ([34]) one has 

{x{t) - f (0)) - e (35) 

upto time r, and ^ t thereafter. 




FIG. 10: Demonstration of the time scaling of Eq. (|35p for 
polymer translocation driven by a field E in 3D. Reproduced 
from Ref. |27l] with permission from Institute of Physics Pub- 
lishing Ltd., UK (original article published on 15 February 
2008). A note about the choice of variables: in order to 
avoid saturation effects [27| the average time (t), required 
to have s monomers translocated for the first time, is plot- 
ted (instead of plotting the average number of translocated 
monomers {s{t)) as a function of t) for N — 400 (average 
over 16, 000 polymer realizations for each field), A'^ = 800 (av- 
erage over 16,000 polymer realizations for each field), and 
N = 1, 200 (5, 000 polymer reahzations for _EA/(fcsT) = 0.05, 
and 7, 500 polymer realizations each for EX/{kBT) — 0.15 and 
EX/{kBT) = 0.25). Here A is the lattice constant. The data 
for N = 800 correspond to real time value, while the data for 
iV = 400 and = 1,200 have been shifted by ^500 units 
along the x-axis for clarity. The solid line has been added for 
a guide to the eye. 

I now demonstrate the time scaling of Eq. ((35)) for 
field-driven translocation of a Rouse polymer, adsorption 
of a Rouse polymer on a solid surface, and translocation 
of a Rouse polymer by means of a pulling force. 
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1. Field-driven polymer translocation and polymer 
adsorption on a solid surface 

With the memory kernel ii{t) for a translocating poly- 
mer shown in Fig. [51 I proceed directly to the simulation 
data to demonstrate Eq. ((35)) for field-driven transloca- 
tion (translocation driven by a field E that only acts on 
the monomer in the pore) and polymer adsorption on a 
solid surface (with energy of adsorption e per monomer), 
both in 3D. These results were originally reported, using 
a Monte Carlo lattice polymer code, in Refs. [23| and 
[28j . where the reader can find the details on the simu- 
lations. The force felt by the monomer in the pore for 
translocation, and by the monomer in contact with the 
adsorbing surface, respectively, are directed perpendic- 
ular to the membrane and the adsorbing surface. Con- 
sequently, their dynamics is described simply by using 
the (scalar) components of the velocities and forces on 
the monomers perpendicular to the membrane and the 
adsorbing surface. 




t 



FIG. 11: (color online) Demonstration of the time scaling of 
Eq. (|35|) for polymer adsorption on a surface. Figure re- 
produced from Ref. with permission from Institute of 
Physics Publishing Ltd., UK (original article published on 12 
May 2009). A note about the choice of variables: in order 
to avoid saturation effects [23|, (instead of plotting the av- 
erage number of adsorbed monomers as a function of time) 
the average number of adsorbed monomers s„ in time t„ is 
plotted, where t„ is the time have n-th monomer adsorbed 
for the first time. Adsorption data with weak adsorption 
energies for A'^ = 1000 (from left to right) £ = 2 (black), 
e = 3 (blue), e — 4 (magenta) and e = 5 (orange): the data 
are progressively separated by a factor 2 along the a;-axis for 
clarity. The original data are shown in the inset in the same 
color scheme. The solid black line corresponds to an exponent 
{l + u)/{l + 2u) ~ 0.73. 

For field-driven translocation in 3D; a polymer of 
length N was threaded fixed halfway through the pore 



and was thermalized with the field E switched on, acting 
only on the monomer in the pore perpendicular to the 
membrane (this makes ip{0) ^ 0). Then translocation 
was started sd, t — 0. It was shown in Ref. \2T\ that 
{(t>{t)) approaches a constant fairly quickly. 

Similarly, for polymer adsorption on a solid surface in 
3D, a polymer of length N was thermalized with one 
end held tethered on the surface. For this case, entropy- 
related stretching of a polymer close to a surface makes 
0(0) ^ 0. Adsorption was started at i — (see Ref. [2^ 
for further details). 

The data for field-driven polymer translocation and 
polymer adsorption in 3D, are reproduced from Refs. [23] 
and [28] in Figs. [10] and [TT] respectively; both confirm 
Eq. ([55)1 . It however needs to be mentioned here that 
for field-driven polymer translocation and polymer ad- 
sorption to a solid surface in 2D, the memory kernel is 
overruled by conservation of energy: for these problems, 
the number of monomers translocated and the number of 
monomers adsorbed respectively scale ^ t^^^^'^\ i.e., the 
translocation and the adsorption time scale as N'^'^ . The 
reader can find more details on this in Ref. UM. 



2. Polymer translocation by a pulling force 

The setup for polymer translocation by a pulling force 
is as follows. A polymer is threaded through a narrow 
pore in a membrane, and a latex bead is attached to one 
of the end monomers. An optical tweezer captures the 
bead, and pulls the polymer perpendicularly away from 
the membrane, facilitating translocation of the polymer 
through the pore [29| . Here 1 consider translocation of a 
self-avoiding Rouse polymer, with the membrane placed 
on the yz-plane, while the bead is pulled with a constant 
force F along the -\-x direction, causing the polymer to 
translocate from left to right. As mentioned earlier, the 
translocation is rate-limited by the events at the pore, so 
the action of the force F alone does not determine the 
translocation dynamics. 

It is clear from this setup that the polymer on the 
right of the membrane would be stretched due to the ac- 
tion of the force, but not on the left of the membrane. 
Therefore, although for this setup the power-law part of 
the memory kernel for the polymer on the left of the 
membrane is still given by ^ one needs 

to establish the memory kernel for the stretched part of 
the polymer on the right of the membrane. Indeed, as 
1 show below, the power-law part of the memory kernel 
for the stretched part of the polymer on the right of the 
membrane is given by t~^^^: it is once again obtained 
from the mean relaxation of a (local) strain in the poly- 
mer at the pore, when the strain is caused by monomer 
insertion into a stretched polymer tethered on a mem- 
brane at the tether point. In order to do so, one does 
need two ingredients: (i) the terminal decay time t of a 
stretched polymer of length iV tethered on a membrane 
scales ~ N'^ |30i], (ii) the shape of stretched polymer is 
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that of a cylinder, where the radius of the cylinder is 
given by that of the Pincus blob ^. The magnitude of 
the pulling force F determines the value of ^, in terms 
of which the spring constant of a stretched polymer of 
length N is given by - ksT/iN^) [J. 

Consider the case when n extra monomers are injected 
at t = 0, at the tether point of a self-avoiding Rouse poly- 
mer, which is tethered on a membrane at one end, and is 
stretched by a force applied at the open end. Given (i-ii) 
above, at time t, counting away from the tether point, all 
the monomers within a backbone distance ^ i^/^ will 
equilibrate to the new situation. However, since the poly- 
mer is stretched, its shape is that of a cylinder, and the 
real space extent of nt monomers is r{nt) ~ nt, but since 
the rest (TV — nt) monomers are not equilibrated to the 
injected monomers at time t, there are (rit + n) monomers 
squeezed in a space that extends only to r{nt). The cor- 
responding compressive force [force = (spring constant) 
X (stretching distance)] from these {nt + ri) monomers, 
felt at the tether point, and hence /i(i), is the given 
by [Sr{nt)]/{nt(,) ~ nl drjut) /dnt]/ (ntO = n/{nt£,), 
which scales ^ [30|. (Once again, this behavior 

lasts only till the terminal time r.) The explicitly eval- 
uated t~'^/'^ power-law behavior of the mean relaxation 
response of a stretched self-avoiding Rouse polymer, teth- 
ered on a membrane, to a local strain caused by monomer 
injection at the tether point is shown in Fig. [T^ The 
data in Fig. [1^ obtained by using a Monte Carlo lat- 
tice polymer model, are reproduced from Ref. [soj, from 
where the reader can find the details of the model. 
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FIG. 12: Mean relaxation response of a stretched self-avoiding 
Rouse polymer, tethered on a membrane, to a local strain 
caused by monomer injection at the tether point at t = 0. 
The data are obtained for = 100 and FX/iksT) = 1, 
using a Monte Carlo lattice polymer model. Here A is the 
lattice constant. The solid line corresponds to the power-law 
gure reproduced from Ref. [30] with permission 
from Elsevier Inc. 

Having put all the above together, one finds that for 
translocation of a self-avoiding Rouse polymer with a 



pulling force, there are two different (power-law) memory 
kernels for the polymer on the two sides of the membrane: 
on the left it is i-(i+'')/(i+2'') and on the right it is i^^/^. 
The system's dynamics is determined by the slower of the 
two: i.e., by t"^/^. Thus, when the polymer is threaded 
fixed halfway through the pore and thermalized with the 
force on one of the end monomers, 0(0) ^ 0, and one 
expects Eq. (|55|) to hold true, with a = 1/2. This is 
verified, using a Monte Carlo lattice polymer model in 
Fig. [131 The data are reproduced from Ref. [sO]- 
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FIG. 13: Mean time required, for FX/{kBT) — 1, to translo- 
cate s monomers for the first time, for s — 5, 10, 15, . . . , A'': 
from left to right TV = 200, iV = 400, iV = 600. The time-axis 
corresponding to N = 400 is the true time, for the A'' = 200 
and A'^ — 600 cases the time axis is divided and multiplied 
by a factor 2 respectively. The angular brackets denote an 
average over 48, 000 polymer realizations. A note about the 
choice of variables: in order to avoid saturation effects [30] the 
average time {t) required to have s monomers translocated for 
the first time is plotted (instead of plotting the average num- 
ber of translocated monomers {s{t)) as a function of t). The 
solid line corresponds to an exponent 1/2. Figure reproduced 
from Ref. (30| with permission from Elsevier Inc. 



V. A CRITIQUE ON THE USE OF 
FRACTIONAL FOKKER-PLANCK EQUATION 
TO DESCRIBE ANOMALOUS DYNAMICS FOR 
POLYMER TRANSLOCATION 

In this short section, 1 return to the issue of non- 
Markovian property of the anomalous dynamics in poly- 
meric systems. 

The fact that in polymeric systems, the motion of a 
tagged monomer is non-Markovian is as such not sur- 
prising: in 3D, a polymer with N monomers resides in 
a phase space of 6iV dimensions (in many-polymeric sys- 
tems such as polymer melts, the dimension of the phase 
space is 6 times the entire number of monomers in the 
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systems) , while the dimension of the phase space associ- 
ated with a tagged monomer is only 6. Thus, if one is to 
follow the dynamics of a tagged monomer, one is essen- 
tially taking a projection of the entire system's dynamics 
on to a few degrees of freedom. While the dynamics of 
the entire system is no doubt Markovian, a projection on 
to a few selected degrees of freedom can indeed render the 
dynamics non-Markovian — this is the very foundation 
of the GLE 

As I have shown in this paper, the dynamics of a tagged 
monomer in a wide range of single-polymeric (phantom 
Rouse, self-avoiding Rouse, self-avoiding Zimm, repta- 
tion, translocation through a narrow pore in a mem- 
brane), as well as many-polymeric (polymer melts) sys- 
tems is robustly formulated by the GLE (IT][2]). It is also 
worth emphasizing that — as apparent from Eq. (|23p — 
in this GLE formulation if the concerned monomer makes 
a move at any time, there is an enhanced chance to undo 
this move in subsequent times; this is where anomalous 
dynamics in polymeric systems stem from. Anyhow, the 
GLE describes trajectories in the phase space. Given 
this, a question that naturally arises is: "how does one 
formulate a probabilistic description of the trajectories in 
phase space for anomalous polymer dynamics?" What I 
am appealing to is that the Fokker-Planck equation is a 
probabilistic formulation of the trajectories described by 
the Langevin equation; so what would be an appropriate 
probabilistic description of the trajectories described by 
Eqs. p1^ . or for that matter Eq. (^5)1 ? 

While for now I will leave this question for future re- 
search here, this section will not be complete without 
mentioning that in a limited context — to describe the 
anomalous dynamics of polymer translocation — frac- 
tional Fokker-Planck equation (fFFE) has been postu- 
lated |!7]. In this equation, an extension of the standard 
Fokker-Flanck equation, anomalous dynamics is a con- 
sequence of introducing power-law waiting times before 
each jump of the concerned particle, as the jump length 
and the waiting times for any jump is obtained from fixed 
probability distributions, independently of their values at 
previous jumps. Given that there is no power-law wait- 
ing time for the movements of the concerned monomer 
(instead if it makes a move at any time, there is an en- 
hanced chance to undo this move in subsequent times; 
this is where anomalous dynamics in polymeric systems 
stem from), not only that the application of fFFE to de- 
scribe polymer translocation is not correct, but also it 
remains a far cry for describing the examples of anoma- 
lous dynamics in polymeric systems that are considered 
in this paper. 



VI. CONCLUSION 

Extending the work of a recent Letter |5|], in this 
pedagogical paper I extensively elaborated on the fact 
that the anomalous dynamics of a tagged monomer in 
a wide range of single-polymeric (phantom Rouse, self- 



avoiding Rouse, self-avoiding Zimm, reptation, translo- 
cation through a narrow pore in a membrane), as well 
as many-polymeric (polymer melts) systems is robustly 
formulated by the GLE. In this GLE formulation the ve- 
locity v(t) of a tagged monomer and the force (/)(t) it 
experiences, are related to each other via 

${t) = - [ dt'fi{t - t')v{t') + g{t). (36) 

JO 

In Eq. (|36p /i(<) is the memory kernel, and the stochastic 
noise term g{t) satisfies the condition that {g(t))o — 0, 
while the fluctuation-dissipation theorem (FDT) {g{t) ■ 
g{t'))o = {${t) ■ (f{t'))ij=o = 3kBTfi{t - t') in 3D. Here 
ks is the Boltzmann constant, T is the temperature, and 
(. . .)o denotes an average over the stochastic noise real- 
izations, including an average over equilibrium configu- 
rations of the polymers at t = 0. Equation (|36p can be 
inverted to write 

v{t) = - [ dt'a{t - t')${t') + h{t), (37) 

with jl{s)a{s) = 1 in the Laplace space, {h{t))o — 0, and 
the corresponding FDT {h(t)-h{t'))o = {v{t)-v{t'))\;j;^^ 
iksT a{t — t'). On the one hand /i(i) is the mean re- 
laxation response of the polymers to local strains, and 
can be derived from the equilibrium statistical physics of 
polymers; and on the other, a{t) characterizes the anoma- 
lous dynamics via the FDT: as the mean-square displace- 
ment of a tagged monomer is obtained by integrating 
{v{t) ■ v{t')) twice in time. An important property of 
the anomalous dynamics that transpires through this ex- 
ercise is that if ^{t) ~ t~°' for some a, then the anomalous 
dynamics is also a. In other words, the anomalous dy- 
namics for polymeric systems are connected to the mean 
relaxation response of the polymers to local strains. The 
anomalous dynamics and the mean relaxation response of 
the polymers to local strains, as shown in Eq. (|23p works 
in the following way: if the concerned monomer makes a 
move at any time, there is an enhanced chance to undo 
this move in subsequent times; this is where anomalous 
dynamics in polymeric systems stem from. 

Further, the characteristics of the drifts caused by 
a (weak) applied field on polymeric systems, too, are 
obtained from the corresponding memory kernels: if 
li{t) ~ t~°' for some a, then the drift exponent is also 
a. This could be thought of as the analog of the Nernst- 
Einstein relation. 

Given that the GLE provides the trajectory description 
in the phase space, I bring to light the non-Markovian 
character of the anomalous dynamics for polymeric sys- 
tems. The fact that the motion of a tagged monomer is 
non-Markovian is as such not surprising: in 3D, a poly- 
mer with N monomers resides in a phase space that has 
QN dimensions (in many-polymeric systems such as poly- 
mer melts, the dimension of the phase space is 6 times 
the entire number of monomers in the systems) , while the 
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dimension of the phase space associated with a tagged 
monomer is only 6. Thus, if one is to follow the dy- 
namics of a tagged monomer, one is essentially taking a 
projection of the entire system's dynamics on to a few 
degrees of freedom. While the dynamics of the entire 
system is no doubt Markovian, such a projection on to 
a few degrees of freedom can indeed render the dynam- 
ics non-Markovian — this is the very foundation of the 
GLE [3ll.[3^. Given this, a question that naturally arises 
is: "how does one formulate a probabilistic description of 
the trajectories in phase space for anomalous polymer 
dynamics?" While for now I leave this question for fu- 
ture research, I note that anomalous dynamics in poly- 
meric systems cannot be captured by, e.g., the fractional 
Fokker-Planck equation (fFPE) which has been recently 
postulated [7] to describe the anomalous dynamics of 
polymer translocation through a narrow pore in a mem- 
brane. In this equation, an extension of the standard 
Fokker-Planck equation, anomalous dynamics is a con- 



sequence of introducing power-law waiting times before 
each jump of the concerned particle, as the jump length 
and the waiting times for any jump is obtained from fixed 
probability distributions, independently of their values at 
previous jumps. Given that there is no power-law wait- 
ing time for the movements of the concerned monomer 
(instead if it makes a move at any time, there is an en- 
hanced chance to undo this move in subsequent times; 
this is where anomalous dynamics in polymeric systems 
stem from), not only that the application of fFPE to de- 
scribe polymer translocation is not correct, but also it 
remains a far cry for describing the examples of anoma- 
lous dynamics in polymeric systems that are considered 
in this paper. 
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Appendix: Derivation of {v{t) ■ v{t')) for phantom polymers [i.e., Eq. (|25p ] 

I start with the velocity autocorrelation function as in Eq. (P5)l 

V{tX)^{v{t)-v{t'))= I dh(3{t-h) f dt2p{t' -t2)[{g{h)-g{t2)) + {fN/2{tl)- fN/2{t2))l (Al) 



A(tl,t2) 



I am interested in the behavior of {v{t) ■ v{t')) in the limit of large {t,t') with finite [t — t'). In this limit I expect 
V{t,t') — >■ V{t — t'). Nevertheless, I take a dual Laplace transform of V{t,t'): one with Laplace variable s for i, and 
the other with Laplace variable s' for t' . 



V{s,s') ^ I dte-'' 



dt' e-'* V{t,t') 



(A2) 



dt'e~''*' / dti(3{t-ti) / dt2P{t' ~t2)[G{ti,t2) + A{ti,t2)]. 
Jo Jo 

Upon having interchanged integration variables, I rewrite V{s, s') as 

y(s,s')= / dte-'* dt'e-^'"' dtiPit-h) dt2 f3{t' ~ t2) [G{h,t2) + A{ti,t2)] 
Jo Jo Jo 

i /"OO pOO POQ 

dhe-'*' dt2e-'''^G{h,t2)+A{ti,t2)] dt e-'^-'-'^^ (3{t ~ h) dt' e'^'^*"*^) /3(t' - ta) 

Jo Jti Jt2 

= ^(s)^(s') / dhe-''^ dt2e-''''[G{h,t2) + A{h,t2)]. (A3) 
Jo Jo 

With G{ti,t2) = SfesT/xd^i - t2\) and A{ti,t2) = &-ikBT5{ti - ia), I choose to represent A{t - t') as 

flimA^os if|i-t'|<A/2 
A{t - t') = 67fcsT <^ 

[ otherwise. 

Thereafter, having recalled that /?(s) = [7 + /i(s)]^^, I obtain 



(A4) 



V{s,s')=ikBTp{s)P{s') 



s' — s s — s' 



s' — s s — s' 



(A5) 
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which implies, assuming t > t', that 

{v{t) ■ v{t')) = SkeTPit - t'). (A6) 
I 
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